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Correlation Effects on Optical Conductivity of FeSi 
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Effects of electron correlation in FeSi are studied in terms of the two-band Hubbard model 
with the density of states obtained from the band calculation. Using the self-consistent second- 
order perturbation theory combined with the local approximation, the correlation effects are 
investigated on the density of states and the optical conductivity spectrum, which are found to 
reproduce the experiments done by Damascelli et al. semiquantitatively. It is also found that 
the peak at the gap edge shifts to lower energy region by correlation effects, as is seen in the 
experiments. 
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FeSi is a well-known material for its anomalous properties that cannot be understood from 
conventional band theories. At low temperatures FeSi is a nonmagnetic semiconductor with a 
gap size of about 50 meV, whereas it has a magnetic moment and shows metallic behaviors at 
room temperatures.!'!) To explain the magnetic properties, Takahashi and Moriya!) applied the 
spin fluctuation theory leading to an idea of thermally induced local moment, which was confirmed 
later by the neutron scattering experiment.!) Recently interest has been renewed by the suggestion!) 
that FeSi can be viewed as a strongly correlated insulator or a Kondo insulator. Kondo insulators 
have been found in the /-electron systems and typical examples are YbBi2,!) CesI^Pts!) and so on. 
They have correlated /-bands and become insulators with small energy gaps at low temperatures. 
The correlation in FeSi may not be so strong as in these materials. It will be, however, of interest 
to reinvestigate the properties of FeSi from a new point of view. 

Low temperature properties of FeSi have been studied intensively from this aspect. Especially, 
Schlesinger et a/.i) measured the optical conductivity and reported that the gap (~50 meV) at low 
temperatures disappears at a temperature (~250 K) lower than the gap size, which they attributed 
to the correlation effect. The optical conductivity of FeSi was reported by several groups. 110.0) 
The gap opening of about 50 meV and the rapid temperature variation were commonly observed. 
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Those optical measurements also provided us with a direct information about the dynamical struc- 
ture of this system at low temperatures such as a precise gap size. Moreover, high-resolution 
photoemission measurement was done by Saitoh et al, who estimated the self-energy of the 
d-electrons and demonstrated a possible existence of the strong renormalization of the electronic 
states at low energies. A recent experiment reported that a sign of the correlation effects is found 
in the tunneling spectroscopy measurements,© where the gap in the density of states (DOS) seems 
to be strongly temperature dependent. It was also reported that the low temperature transport 
measurements for Fei_ x Co x Si and FeSii_ x .AlJll) suggest strong mass enhancement of carriers. 

In these contexts several theoretical attempts have been made to take the correlation effects into 
account. Among these attempts Fu and Doniach0) proposed a two-band Hubbard model based 
on their band calculation. Their calculation revealed the importance of the correlation effects 
on physical quantities but there seems to be some errors about the treatment of the self-energies. 
Using this model, we© also calculated physical quantities correctly, and showed that the tem- 
perature variations of the magnetic susceptibility and the optical conductivity are consistent with 
the experiments. From a viewpoint of the strongly correlated systems, Riseboroug 

hi) 

calculated 

the angle-resolved-photoemission spectrum of the Anderson lattice model in terms of the 1/N 
expansion. 

On the other hand several band calculations have been done for FeSi. Each of them 

succeeded in giving the insulating ground state and showed that the minimum gap is given by the 
indirect one, but it is close to a direct one. It has been also demonstrated that the states near the 
Fermi level mainly consist of the d-electrons. However, somewhat different gap sizes (0.03-0.2 eV) 
were reported. There are some attempts to explain the physical properties of FeSi by extending the 
band calculations, e.g., Jariborgil) investigated the thermal disorder effect and Anisimov et al^ 
applied the LDA+C7 scheme. Fu et al& calculated the optical conductivity without the correlation 
effect using their band calculation. As a result, the necessity of the correlation effects was underlined 
in order to reproduce the observed temperature variation. Ohta et a3 also calculated the optical 
conductivity in the joint-DOS form from their band calculation and showed that the gap is filled 
with increasing temperature. In their calculation, the Drude contribution is estimated from the 
experimental DC resistivity and is added to the joint-DOS term separately. However, the flat part 
of the optical conductivity spectrum within the gap is not reproduced. 

Recently, Yamada et a/.© calculated the DOS within the LMTO-ASA method carefully, and 
obtained the gap size close to the experimental one. Therefore, in the present Letter, we study the 
two-band Hubbard model using the Fe part of the DOS obtained from this band calculation and 
calculate the temperature variation of the optical conductivity to investigate the correlation effects. 
Moreover, we compare the calculated spectrum of the optical conductivity with the experimental 
data and estimate the gap size (indirect one) in the DOS of the real system. 
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We start from the following two-band Hubbard Hamiltonian: 



i 



+ U 2 ^( n iiT n i2| + n i 2^n ill ) 



i 



+ U 3 ^2(ni^n i2 ^ + n i2 in in ) 



i 



(i) 



where the (cj ao -) creates (destroys) an electron on site i in band a = 1,2 with spin a. If the 
two bands are degenerate, this Hamiltonian is rotationally invariant in spin and real spaces when 
one chooses U 2 = U — J and U3 = U — 2jjH3 where U and J denotes the Coulomb repulsion 
and the exchange interaction, respectively. We adopt this relation so as to reduce the number of 
parameters. 

The initial DOS obtained from the band calculations' is displayed in Fig. 1 (see the solid line 
for T = 0), where we widen the energy gap by 16 % so as to reproduce the gap edge of the optical 
conductivity at 4 K of the experiment. There is an ambiguity to determine the gap size E g . If 
we extrapolate the steepest part of the DOS at the both sides of the gap neglecting the tails, we 
obtain E g = 75 meV, whereas if we regard the gap as the region inside the tails of the gap edge, we 
obtain 60 meV. The Fermi level {Ep = 0) is placed at the center of the gap at T = 0, and we call 
the upper and the lower part of the DOS as the band 1 and 2, respectively. We introduce a cut off 
for each band so as to include one state per spin in each band. Then the band width for the band 
1 is about 0.56 eV and about 0.63 eV for the band 2. Moreover, although the DOS is asymmetric, 
we assume that the chemical potential is temperature independent. This means that the electron 
number of this system varies with temperature. This assumption seems to be reasonable in the 
real system, where electrons and holes are supplied from other bands. 

Because the present interests are on the correlation effects in the low energy and low temperature 
region of this model, we adopt the self-consistent second-order perturbation theory (SCSOPT) 
combined with the local approximation.® In SCSOPT, the self-energies (uS) and (cj) are 
calculated as 
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where 



and 
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Here, N is the number of sites, /(e) the Fermi function and p^{e) the DOS of band a for the 
non-interacting case. We solve the set of equations (2)-(4) self-consistently. To make numerical 
calculation easy, we take 5 finite (5 = 10 -4 ) in eq. (|3|) and convert these equations with the 
transformations© 

/oo 
dee~^pl(e)f{e), 
-00 

B° a (r) = / dee-^pl(e)f{-e). (5) 



In general, the current operator is expressed as 

i = e££< m ' c Lk<Wk, (6) 

cr,k mm' 

where m denotes the band index. To derive the expression for the optical conductivity, we assume 
the intra- and interband contribution to be equal (v^ m ' = Using the linear response theory, we 
assume that the momentum conservation is violated in real systems by defects and phonon-assisted 
transitions and replace ^k^k'^kk' with v 2 /N in the current-current correlation function. Then, we 
finally obtain a joint-DOS-like form for the optical conductivity: 



0{u,T) = — - — 2^ / du 



ft a J -oo oj 

x [pU") + pliy)] [pi iy + + fiiy + u)]. (7) 

In the numerical calculation we set (ev) 2 jfi = 1 for simplicity. 

Figure 1 shows the temperature dependence of the quasiparticle DOS for U = 0.5 eV and J = 
0.35U. These values are chosen so as to reproduce the shape and the temperature dependence of 
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Fig. 1. The initial DOS obtained from the band calculation(Ref. 26) at T = 0. At finite T, the DOS is strongly 
temperature dependent due to the correlation effects. 



the optical conductivity spectrum. At K, the band 1 is filled and the band 2 is empty, so that 
the correlation effect is absent except the Hartree-Fock contribution. At finite T, however, the 
correlation is introduced through the thermally excited electrons and holes. As a result, the DOS 
becomes strongly temperature dependent and is smeared especially near the gap. The gap at K 
is almost filled up at the temperature of the order of its size and it is also found that the peaks of 
the gap edges move toward the center of the gap. 

In Fig. 2(a), the temperature variation of the optical conductivity calculated from eq. (7) is shown. 
Since the intraband (Drude) contribution does not exist at K, only the interband contribution 
survives and seems to reproduce the shape of the spectrum of the experiment at 4 K. The slightly 
rounded gap edge in the calculation is due to the tails at the gap edge in the DOS of the band 
calculation.© If we extrapolate the steepest parts of the edge, we obtain the gap size as about 75 
meV, which is consistent with the gap size of the DOS. Since these features reflect the shape of 
the DOS directly (see eq.(7)), the band calculation using LMTO-ASA by Yamada et aZ0 seems 
to give good results about the whole structure of the DOS at T = except for a slightly smaller 
gap size. 

At finite temperatures, the experiment done by Damascelli et al& (Fig. 2(b)) shows not only 
that the gap is almost filled up at 300 K, but that the increase of the weight in the gap from 100 to 
300 K seems to be more rapid than that from 4 to 100 K. One can also recognize that the peak at 
the gap edge shifts to the lower frequency region as the temperature rises. In our calculation (Fig. 
2(a)), the gap is almost filled up at 300 K as well as the rapid increase from 150 to 300 K, but a little 
slower compared to the experiment, which may be mainly because the absence of the correlation 
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Fig. 2. (a) The temperature dependence of the optical conductivity calculated with the eq. (7). (b)The experimental 
data from Ref. 11. The peaks due to phonons observed in the gap are omitted. 



effect at K in the present model. Reflecting the correlation effects the peak at the gap edge shifts 
to lower frequency region, as is seen in the experiment. In our calculation, however, there are dips 
between the Drude and the interband contributions in contrast to the experiment. This may be 
caused by the assumption in the previous formulation, where we neglected the k-dependence of 
the current operator and took the intraband (Drude) and interband contribution to be equal in 
deriving eq. (7). 

For comparison, the optical conductivity obtained from the Hartree-Fock approximation or a 
rigid band model is shown in Fig. 3. Since the gap is filled only with the Drude contribution, the 
temperature variation is monotonous and the gap is not efficiently filled at a temperature of the 
order of the gap size as pointed out by Fu et al& The peak at the gap edge does not shift. These 
results are consistent with the previous studiesji3!) but are in disagreement with the experiments. 

In summary, we have calculated the correlation effects using the two-band Hubbard model and 
the DOS obtained from the band calculation, and compared the calculated optical conductivity 
with the experiment. Although we used an essential but simplified model and applied the local 
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Fig. 3. Calculated optical conductivity within the Hartree-Fock approximation. 



approximations in deriving eq. (7), the calculated optical conductivity reproduced the experiment 
rather well. Thus we determined the gap size (indirect one) as 75 meV when the tail is neglected (60 
meV within the tail). Furthermore, we confirmed that the correlation effect is essential to explain 
the rapid temperature dependence of the optical conductivity, and gave a direct explanation for 
the shift of the peak at the gap edge. The latter is another possible sign of the correlation effects. 
Similar effect is also seen in the optical measurement of YbBi2,© but the shift is observed on the 
IR peak instead of the peak at the minimum gap edge. 

In the present model the higher order terms of the electron correlation effect are quenched at 
T = 0. The renomalization effects at T = 0, however, may be also important as in the model by 
Fu and Doniach® The optical conductivity of this model is strongly temperature dependent as 
is seen in our previous calculation. It is also important to take the spin fluctuations!' U into 
account at higher temperatures or for magnetic properties. Furthermore, to incorporate the details 
of the band dispersions into the methods of the strongly correlated electron systems is our present 
interest. 
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Figure Captions 

Figure 1: The initial DOS obtained from the band calculation(Ref. 26) at T = 0. At finite T, the 
DOS is strongly temperature dependent due to the correlation effects. 

Figure 2: (a) The temperature dependence of the optical conductivity calculated with the eq. (7). 
(b)The experimental data from Ref. 11. The peaks due to phonons observed in the gap are omitted. 
Figure 3: Calculated optical conductivity within the Hartree-Fock approximation. 
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